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Abstract 

A method of exciting coherent spatial modes of Bose-condensed trapped 
atoms is considered. The method is based on the resonance modulation of the 
trapping potential. The population dynamics of coherent modes is analysed. 
The method makes it possible to create mixtures of different spatial modes in 
arbitrary proportions, including the formation of pure excited coherent modes. 
Novel critical effects in the population dynamics are found. 
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1 Introduction 



Bose atoms in magnetic traps can be cooled down to ultralow temperatures where 
Bose-Einstein condensation takes place. There exists numerous literature, both 
theoretical and experimental, devoted to this subject (see reviews [1,2]). 

The Bose condensation occurs when the thermal wavelength of atoms, A = 

27r fi^ /rriokBT, becomes larger than the mean interparticle distance, a. At the 
same time, the effective radius of an atom, ao, has to be much smaller than the in- 
teratomic distance in order that strong hard-core repulsion would not disturb much 
the motion of colliding particles. When the thermal wavelength is such that 

A » a » ao , 

atoms are mutually correlated and a coherent state develops. 

In the system under equilibrium, condensed atoms are in the ground state. Re- 
cently [3] , the idea was advanced of the possibihty to create non-ground-state con- 
densates corresponding to excited coherent modes of trapped atoms. This can be 
done by imposing, in addition to the trapping potential, a time-dependent field 
whose oscillation frequency is adjusted to be in resonance with the transition fre- 
quency between the ground-state level and a chosen excited level. In this report, we 
present the results of investigation of the population dynamics in a coherent system 
of trapped atoms, driven by a resonance field. As it turned out, the dynamics of 
such a system is quite nontrivial and is rather different from the population dynam- 
ics of two-level optical systems. Thus, unusual effects, that can be called critical, 
are found when a small variation of system parameters results in a drastic change 
of the population dynamics. 



2 Population Dynamics 

Consider the Bose gas of neutral atoms at ultralow temperatures, when all atoms 
are in a condensed state. The wave function of coherent atoms is obtained from the 
nonlinear Schrodinger equation 

dip 



dt 



H{ip) + Vres (1) 



in which H{(p) is a nonlinear Hamiltonian and V^es is a resonance field. The gas is 
assumed to be diluted, and the interatomic interactions are shape independent and 
can be modelled by the contact potential 

$(r ) = (5(r ) , 

mo 

where uiq is atomic mass and is the scattering length. Then the nonlinear Hamil- 
tonian is 

H{ip)= -^^' + U{f)+m\ (2) 
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where C/ is a trapping potential and 

A = Anh'' ^ N , 
mo 

represents the effective interaction for N confined particles. The resonance field can 
be taken in the form 

Vres = V{f) COSLUt . (3) 

Let us note that the nonlinear Schrodinger equation is an exact equation of 
coherent states [4] . Contrary to this, if (/? is treated as the order parameter associated 
with the condensate, then Eq. (1) is an approximate equation corresponding to the 
mean-field approach at zero temperature. Such an approximate equation is often 
called the Gross-Ginzburg-Pitacvskii equation [5-7]. 

The stationary solutions of Eq. (1), if the time-dependent field is absent, are 
ipn^w{~iEnt/h) with ipn and En defined by the eigenvalue problem 

H{(Pn)<Pn = EnVn ■ (4) 

The eigenfunctions labelled by a multi-index n, are stationary coherent modes. 
The transition frequencies cUmn between two energy levels are given by the difference 

^t^mn = E^i — En ■ (5) 

Assuming that at the initial time t = 0, all atoms are in the ground-state 
coherent mode, we have the initial condition 

^{f,0) = Mr). (6) 

And let then the resonant field (3) be switched on, with a frequency cu being in 
resonance, or almost in resonance, with the transition frequency 

Ej - Eq 

^jo= (7) 

between the ground state and a chosen energy level j. The corresponding quasires- 
onance condition is 

< 1 , Au; = u;-u;jo. (8) 

The solution of the time-dependent equation (1) can be presented as an expan- 
sion 

ipif,t) = Yl CnitMf) exp (^-^ Ent^ (9) 

over the coherent modes (pn{r)- Combining Eqs. (1) and (9), we find a system of 
equations for the coefficients c„(t). This system can be simplified in the resonance 
approximation giving 



— ^ -I arijCo - 2 ' 



^ = -i an,c, - \ P*e-'^^' c, , (10) 

where 

n,{t) = |Q(t)|2 (11) 

is the population of the level i, and a new notation is introduced for the interaction 
amplitude 

a^j J \Mr)\'Mr)\Uf (12) 
as well as for the transition amplitude 

I I ^oir)V{f)^,{f)df. (13) 

Because at the initial time all atoms are in the ground state, the initial conditions 
to Eqs. (10) are 

Co(0) = l, c,(0)=0. (14) 

Since the quantities Cq and Cj are complex, Eqs. (10) are to be complemented by the 
equations for the complex conjugate Cq and c* or by the equations for the populations 
(11). The latter equations are 

^= Im(/5e^^-cSc,) , 
^ = Im {P*e-'^-' c*Co) , (15) 

with the initial conditions 

no(0) = 1 , n,(0) = 1 . (16) 

In addition, the normalization condition 

no{t) + nj{t) ^ 1 (17) 

holds true. The derivation of Eqs. (10)-(17) was expounded in detail in Ref. [3]. 

The system of nonlinear differential equations (10) and (15) can be solved ana- 
lytically by means of the averaging technique [8] yielding for the populations 

^0 = 1 - ^ sm y , rij^ sm — , (18) 

in which Q is an effective frequency given by the equation 

^2 ^ _ _ Aa;]2 ^ |^|2 _ ^^9^ 

The quantity Q is an equivalent for the Rabi frequency, although one has to keep 
in mind that it is actually not a fixed frequency but a function of time defined by 
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Eqs. (18) and (19). The solutions (18), for < a, give a qualitative picture of 
complicated nonlinear oscillations of the populations (11). In the limit \/3\ ^ a, 
these solutions become asymptotically exact. 

Equations (15) show that if in some moment of time the pumping resonance field 
(3) is switched off, then after this the populations Uq and rij remain constant at their 
instantaneous values. This suggests the way of creating mixtures of coherent modes 
with different spatial configurations, and even pure excited coherent modes. Such 
excited coherent modes will not, of course, last for ever. Their lifetime is defined by 
the total lost rate jj caused by atomic collisions. The same, actually, concerns the 
ground-state condensate whose lifetime is defined by the corresponding lost rate 70. 
The lost rates due to binary collisions can be presented as 

13 = ^33^"^^] j \^j{r)\'^df+\jQN'^njnQ j \Lpj{f)\'^\^Q{r)\^dr, 

^o = XooN^nl J \ipo{f)\Uf+XojN^nonj J |(^o(r )|>,(f )|2 , (20) 

where Xij are the corresponding relaxation coefficients. The lost rates due to ternary 
collisions can be presented in the similar way. The ternary loss rates are usually 
much smaller than the binary ones [9]. 

The approximate solutions (18) and (19), obtained by the averaging method, 
give us a general understanding of the excitation procedure. However, Eqs. (18) 
and (19) define the level populations not explicitly but through a connected system 
of equations. In order to study the dynamics of the populations (11) explicitly, and 
also not to be limited by the averaging-technique approximation, we return back to 
the evolution equations (10), which we solve numerically. For numerical analysis, 
it is convenient to measure time in units of and to introduce dimensionless 
quantities 

a a 

The results of numerical calculations for several values of parameters are presented 
in Figs. 1 to 6. 

It turned out, that the dynamics of the populations is rather nontrivial and 
exhibits a kind of critical effects when an anomalous coherent oscillation of the 
state populations suddenly appears. The latter occur on the critical line defined 
approximately as be = 0.5 — 6. Outside the critical line, the populations oscillate 
according to the sine-squared law, as in Fig.l. Approaching the critical line, for 
instance by changing the detuning, the period of oscillations suddenly increases 
by a jump, with the top of Uj and, respectively, the bottom of no, becoming fiat. 
Thus, the slight change of the detuning, from 6 = 0.11 in Fig.l to 6 = 0.11017478 
in Fig.2, results in an abrupt period doubling. And the following small shift of the 
detuning to6 — 0.11017479 results again in the period being approximately doubled, 
with the qualitative change of the time-dependence: Each second downward cusp 
between two adjacent oscillations of rij overturns up becoming an upward cusp. 
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which yields to an effective period doubling, as is seen from comparing Figs. 2 and 
3. The further increase of the detuning drives the system away from the critical line 
resulting in the decrease of the oscillation period and in the oscillation behaviour 
again more resembling the sine-squared shape, as in Figs. 4 and 5. The overall 
picture is qualitatively the same, when we cross the critical line at other values of 
the parameters b and 6. Thus, Fig. 6 shows the time behaviour of the populations 
at 6 = 0.499 and 6 = 0.001001002 being on the critical line. Figure 6 is similar to 
Fig. 2, except that the oscillation amplitudes are different, being equal to different 
values of b. 

In practice, the value of b, which is the dimensionless expression for the transition 
amplitude (13), depends on the resonant field (3). The spatial part of the latter can 
be chosen to have different forms. For example, in the consideration of a mixture 
of two Bosc condensates of ^^Rb atoms in two internal hyperfine states, en effective 
potential forcing the excitation of a first antisymmetric mode was linear [10], due to 
the spatial separation of the condensate components. In our scheme, the potential 
V{f) can be arbitrary. 

3 Conclusion 

We presented the analysis of temporal behaviour of spatial coherent modes excited 
by a resonant field realizing the oscillatory modulation of the trapping potential. 
Such a procedure of using the resonant pumping field suggests the way of creating 
mixtures of coherent spatial modes in arbitrary proportions, including the formation 
of pure excited coherent modes. 

There exists the critical line connecting the values of the transition amplitude 
and detuning at which the dynamics of fractional populations suddenly changes. 
The qualitative and sharp change in the population dynamics reminds critical phe- 
nomena occurring in equilibrium systems. It is possible to construct an effective 
stationary system describing the averaged behaviour of populations satisfying the 
original evolution equations. The effective averaged system displays critical be- 
haviour at the critical line which corresponds to that observed for the nonequilib- 
rium system. Effective critical indices can also be defined. The detailed study of the 
critical dynamics is presently under investigation and will be published elsewhere. 
The possibility of exciting various coherent modes of Bose atoms can be important 
in the context of realizing different spatial modes of atom lasers [11-15]. 
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Figure captions 



Fig.l. The populations of the excited coherent mode (sohd hne) and of the 
ground-state mode (dashed hne) as functions of time for b — 0.4 and 6 — 0.11. 

Fig.2. The critical dynamics of the populations for b = 0.4 and 6 = 0.11017478. 

Fig.3. Overturning of the downward cusps of rij upward, with the period of 
oscillations being approximately doubled, occuring ai b — 0.4 and S — 0.11017479. 

Fig. 4. The time dependence of the populations outside the critical line, with 
b = 0.4 and 6 = 0.12. 

Fig.5. The fractional population oscillations for b — 0.4 and 6 — 0.2. 

Fig. 6. The critical population dynamics for b = 0.499 and 6 = 0.001001002. As 
in all previous figures, solid line corresponds to rij and the dashed one to no. 
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